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Abstract 

The Schrodinger-Coulomb Sturmian problem in R , A^ ^ 2, is considered in the momen- 
tum representation. An integral formula for the Gegenbauer polynomials, found recently by 
Cohl arXiv: 1105.2735 , is used to separate out angular variables and reduce an integral Stur- 
mian eigenvalue equation in R^ to a Fredholm one on R+. A kernel of the latter equation 
contains the Legendre function of the second kind. A symmetric Poisson-type series expansion 
' of that function into products of the Gegenbauer polynomials, established by Ossicini [Boll. 

Un. Mat. Ital. 7 (1952) 315], is then used to determine the Schrodinger-Coulomb Sturmian 
eigenvalues and associated momentum-space eigenfunctions. Finally, a relationship existing 
between solutions to the Sturmian problem and solutions to a (physically more interesting) 
energy eigenvalue problem is exploited to find the Schrodinger-Coulomb bound-state energy 



(N 
> 



levels in R , together with explicit representations of the associated normalized momentum- 
\£) ' space Schrodinger-Coulomb Hamiltonian eigenfunctions. 

KEY WORDS: Coulomb problem; momentum representation; integral equations; Gegen- 
bauer polynomials; Legendre functions 
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1 Introduction 

A closed-form expression for the bound-state hydrogenic wave functions in the momentum repre- 
sentation in R 3 was first obtained in 1929 by Podolsky and Pauling p], who succeeded to carry 
out a Fourier transform of the position-space Coulomb eigenfunctions, found three years earlier 
by Schrodinger 2 . Guided by the results of Ref. [T], in 1932 Hylleraas 3 derived a second-order 
partial differential equation obeyed by the momentum-space Coulomb wave functions. In 1933, 
Elsasser [1] attacked the problem in the spirit of Ref. (T], but with the Fourier transform done 
with the aid of the complex contour-integration technique. A completely different approach to the 
subject was presented in 1935 by Fock [5]. He showed that the momentum-space Schrodinger- 
Coulomb wave equation in R 3 , known to be an integral one, might be transformed into a form 
which appeared to be identical with an equation satisfied on the unit hypersphere § 3 by the hy- 
perspherical harmonics. This fact led Fock to the conclusion that for bound states the symmetry 
group of the quantum-mechanical Coulomb problem in R 3 is 0(4). That thread was pursued fur- 
ther a year later by Bargmann [6] . Since then, the problem was revisited and discussed in various 
aspects by a number of authors. Space limitations allow us to mention here only several represen- 
tative papers [7H20j. as well as relevant reviews and monographs [2 1- 28] , where references to other 
pertinent works may be found in abundance. 

* Email: radek@mif.pg.gda.pl 
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In this paper, we present a method which allows one to determine bound-state energy levels and 
associated Hamiltonian eigenfunctions of the Schrodinger-Coulomb problem in the iV-dimensional 
(N ^ 2) Euclidean momentum space, making no use of the 0(N + 1) dynamical symmetry specific 
to the problem nor of the Erdelyi's generalization [29] (cf. also Ref. [30l Sec. 11.4]) of the Funk- 
Hccke theorem [311132] regarding the integral equation obeyed by the hyperspherical harmonics 
(which constituted the basis of the approach presented by Levy [7] in the case of N = 3). We 
proceed along the following steps. First, we replace the integral energy eigenproblem by a related 
Sturmian integral eigenvalue equation, in which energy is fixed and the role of the eigenparame- 
ter is played by the Coulomb potential strength. Next, we reduce the aforementioned Sturmian 
integral equation in to a Fredholm one on R + , exploiting a recent result from the theory of 
Gegenbauer polynomials due to Cohl [33] . A kernel in the latter one-dimensional equation contains 
the Legendre function of the second kind. We solve that integral equation using a Poisson-type 
expansion of the aforementioned Legendre function in Gegenbauer polynomials, derived long ago 
by Ossicini [34] . but now buried deeply in the mathematical literature and little known. Finally, 
we use a relationship existing between solutions to the Sturmian problem and to the energy eigen- 
value problem, and determine the Schrodinger-Coulomb bound-state energy levels in M. N , together 
with explicit representations of the associated normalized momentum-space Schrodinger-Coulomb 
Hamiltonian eigenfunctions. 

2 Integral equation for the momentum-space Schrodinger- 
Coulomb wave functions 

In the position representation, the time-independent Schrodinger equation for a nonrelativistic 
spinless particle moving in the iV-dimensional (N 2) Euclidean space M. N , in the field of a force 
derivable from the scalar local potential V(r), is 

h 2 

- —V 2 t>(r) + V(r)^(r) = E$(r). (2.1) 
2m 

The counterpart equation in the momentum representation is 

(s - E ) *<"> + idjm I 'V - #M = < 2 ' 2 > 

where $(p) and U(p) are the Fourier transforms of ^>(r) and V(r), respectively. In the case of the 
attractive Coulomb potential 

V ^ = -7T^Z ( Z> °)' ( 2 - 3 ) 
(47re )r 

which concerns us here, U (p) may be shown (cf. Appendix) to be 



Ze 2 {2h) N ' 2 - l T (Ofl) 1 
(47re ) y/n p T 



so that Eq. (|2.2[) becomes 



t- - e) *( V ) Ze " r( ^ } [ d N v< (2 5) 

2m E ) nP) (Ane ) 7 R « d P \p> p\^ ' (2 ' 5) 

We aim at finding solutions to the above equation which are associated with, yet unknown, negative 
(bound-state) energy eigenvalues E. 



2 



3 Solution of the integral equation for the momentum-space 
Schrodinger— Coulomb Sturmian functions 

3.1 Preliminaries 

Instead of attacking Eq. (|2.5[) directly, it is easier to consider first the associated Sturmian eigen- 
problem 

— - Ej Z(E,p) - A (47r£o) 2n(N+1)/2h d P lp ,_ p \ N - 1 > (3-1) 

in which A = X(E) is an eigenvalue, whereas E < is a fixed parameter. It is evident that once 
the A-spectrum ({A a }) for Eq. p.l[) is determined, the sought negative energy eigenvalues to the 
original problem (|2.5[) are roots to the algebraic equations 

Xa(E) = 1 (E< 0). (3.2) 

If E a p < is a particular solution to Eq. (|3.2[) (the second subscript at E serves to distinguish 
between different roots to Eq. (|3.2[) . if there are more such roots than one), the corresponding 
negative-energy momentum-space Coulomb wave function $ Q ^ 7 (p) is 

$<*/3 7 (p) = AZ aj (E aP ,p), (3.3) 

A being an arbitrary non-zero constant (the second subscript at S appears if the eigenvalue A a is 
degenerate) . 

Before proceeding further, we shall rewrite Eq. (13. 1[) in a slightly more compact form. To this 
end, we define 

q = y/-2mE, q B = — (3.4) 
a B 

(here a B = (47reo)fi. 2 /me 2 is the Bohr radius). Use of Eq. (|3.4j) transforms Eq. (|3.1j) into 



^ + f)H(E,p) = X ZqB J^ dV W \ p \ N -, - (3-5) 

3.2 Reduction of Eq. ( 13. 5 j) to the radial form 

We try whether Eq. (|3.5p has solutions in the form 

X lv (E,p) = F l (E,p)Y l ^- 1 \n p ). (3.6) 
Here Y^'^ (n p ), with I e N and r\ € {1, 2, . . . , d, (JV-1) }, where 



are the hyperspherical harmonics |25U26j . that are orthonormal on the unit hypersphere 
the sense of 

d^ 1 ^ y^- X> {n v )Y^r X \u p ) = 5 lv 8 m , (3.8) 

We insert Eq. (|3.6[) into Eq. (13.5[) . multiply the resulting equation by Y i ^ V ~ 1 ^*(n p ), and then sum 
over r). Using the addition theorem [26j Eqs. (3.80), (3.86) and (3.40)] 

d ' N ~ 1} ,(Jf-l) r (JV/2-l) (n ,x 

V y M 'fn ^fn')-^ ' — p) (39) 

2^ Y iv y n p) Y i v ^ n p> - s cr (JV/a_1) (i) > ^-v) 
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where Cn"\£) is the Gegenbauer polynomial and S/v-i is the surface area of S w 1 , we obtain 
(p 2 + q 2 )F l (E,p) = XZq B 1 2 I n / dp'p'^-^F^E^') 

d N -i n > p — L -i — [np y N (3.io) 

It is evident that the angular integral appearing in Eq. (|3.10[) may be expressed as 

d ^-i n ; — h (np y N = s N - 2 / d£(i-n (JV ~ 3)/2 — 3 «) . 

(3.11) 

However, it has been recently shown by Cohl [33J that 

1 r .c ( 1 ~^ 2 )" 1/2g " JO _ i7 r (a- K )oa+l/2 r (q)/ 1 N r ( a + |) / 2 ■. \ (q-«)/2fl«-q / \ 

./ 4 (z _ C)K+ l/2 -° 2 °« W r ( K +|)^ ' U n+a-l/2^ 

(Rea > —\i k € C, n e N, z eC\(-oo,l]), (3.12) 

where 0£(z) is the associated Legendre function of the second kind. It follows from Eq. (|3.12l) 
that 

\ (1 " ( fr alvf (0 = 2 Q+1 / 2 ^)(i)a l+Q - 1/2 w 

(Rea > -i, n e N, z G C\(-oo,l]) (3.13) 
(0,/(z) = is the Legendre function of the second kind), and consequently 



d iv ~v 



- - rip ■ n p J 

(3.14) 



2pp> p J 

Plugging Eq. (|3.14p into Eq. (|3.10p and making use of the relation 



(TV -l)^- 1 )/ 2 

r(^±i) 



^-2 - ' ( 3 - 15 ) 



we arrive at the following integral equation for the radial function Fi {E, p) : 

(P 2 + q 2 )F l (E,p) = X^^[dp'p'^ 2 Q l+(N ^ 2 (^f) Fl(E, P ') (3.16) 

(for N — 3 and A = 1, Eq. (13.16[) coincides, as it should, with Eq. (8.7) in Ref. [35). After some 
straightforward rearrangements, this may be further cast into the form 



p 2 +p' 2 



X 2 -^ dp' -==y-^ ^[,/p7^ p '(N-X)/2 Fl{EM (3<17) 



yjp 2 + q 2 \Jp' 2 + q 2 
the kernel in the above equation being real and manifestly symmetric. 
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3.3 Intermezzo on the spectral theory for Predholm integral equations 

Before we attack Eq. (|3.17l) . we have to recall some very basic facts from the spectral theory for 
integral equations. If the kernel M(p,p') in the homogeneous Fredholm equation of the second 
kind 

poo 

f(p)=X dp' M(p >P ')f(p') (3.18) 
Jo 

is real and symmetric, and if its A-spectrum (known to be real, in view of the aforementioned 
assumptions about the kernel) is purely discrete, the following spectral expansion of M(p,p') 
holds: 

M(p,p') = ^2 fk(p \ fk(p ' ) , (3.19) 
k Xk 

provided the eigenfunctions {fkip)} are orthonormal in the sense of 

dpf k (p)fk'(p) = S kkl . (3.20) 



This implies that if we cope with some particular equation of the form (|3.18[) . and if, in what- 
ever manner and from whichever premises, we infer that the kernel M(p,p') possesses the series 
representation 

M(p,p') =J2 £ k9k(p)gkip') (e* = ±1), (3.21) 

k 

the functions {gk(p)} being orthogonal, 



dpg k (p)gk'{p) = \\gk\\ Skk>, (3.22) 

we may immediately deduce that the eigenvalues for Eq. (|3.18[) arc 

Afc =e fe ||g fc |r 2 , (3.23) 
whereas the associated normalized eigenfunctions are 

fkip) = fr^TT = V s k^k9kip) 



Wk\ 

(observe that e^ 1 — £&). 

3.4 Solution of the radial integral equation (I3.17j) 

We return to Eq. (|3.17j) . It is evident that it falls into the category (|3.18j) . with 



(3.24) 



flip) = Vp T +^P {N ~ 1)/2 F 1 (E,p) (3.25) 
(here the subscript I at / corresponds to the one at F) and 







97n„ ~ l + ( N -3)/2 , „ , 

Ml (p,p>) = ^ 2 \ t ] PP 2 / - ( 3 - 26 ) 

7T VP + 9 V -P + 9 

An expansion of the kernel Q3.26P in the symmetric form (|3.21[1 may be derived from the following 
Poisson-type identity: 

qJ ) = 2 ^ir(„+i)}^i-ef +1)/2 ii-e 2 f +1)/2 



p2 _|_ pi 2 



2hy/(i -e 2 )(i-e' 2 ) 



7l r —0 

(3.27) 
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discovered several decades ago by Ossicini [Ml Eq. (11)] (there is a misprint in the formula provided 
originally by Ossicini — on its left-hand side [r(A)] should be replaced by [r(A)] 2 ; moreover, one 
should be aware that in Ref. [361 Eq. (47.6.11)] the formula in question was reprinted with several 
errors!). For h = 1, Eq. (|3.27[) becomes 

qJ m =] = 2^[r { , + i)r(i-ef +1)/2 (i-rf +1)/2 



V — C (iy+1) (£)C^ +1) (£') 

^T(n T +2v + 2) n * lCj n ' {Z) 

(3.28) 



from which, with the substitutions 



^2 — p2 ^ ^2 — p/ 2 

£ = — 5~" i 2 ' £ = ~ 2~ i /¥' (3.29) 



we infer that 



2pp' ) \ q 2 + p 2 J \q 2 + p' 

x E r< It + 2^ (4x4) (3-30) 

r(rt r + 2v + 2) \q +P J \q z +p' 2 J 

(the particular reason for which we have denoted the summation index as n r will become clear 
shortly). Combining Eqs. (|3.26[) and (|3.30p . we see that the kernel in Eq. (13.17)) may be written as 

oo 

Mi{p,p') = 9nAp)9nAp') (=*• Vn r G N : e ril i = 1), (3.31) 

n T =0 

with 



In the next step, we evaluate the integral J °° dp 9n r i{p)gn' T i(p), in order to verify whether the 
functions (13.32[) form an orthogonal set [as they should, cf. Eq. (|3 . 22[) ] . Exploiting the known [371 
Eq. (7.313)] integral formula 

« (i - er-^c^mo - 2 ,-,;s:;;^ (a)P ^ <^ > -i, ^ ■>>. 

(3.33) 

we hnd that 

/ dpg ni i{p)g n[ i(p) = : WZT~ S n r n' t , (3-34) 

Jo n r + I H — q 

i.e., the functions (|3.32[) are indeed orthogonal, the norm of g ni i(p) being 



Z q B 
n r + 1 + 



\9n t l\\ = xl , , ; N _ ± — ■ (3.35) 



Hence, applying Eqs. (|3 . 23[) and (13.24[) we infer that the spectrum of eigenvalues in Eq. f|3 . 1 7(1 is 

given by 

n + 1 + E=l n 

Ki= 7 2 ~ (nreN), (3.36) 
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whereas the associated orthonormal eigcnfunctions are 



f nr i(E,p) = Til 



N-l\ qn T \(n x + l 



N-l - 
2 ' 



n(n T + 2l + N -2)! 



{Aqp y + (N-i)/2 _ 1)/2) 

( q 2 +p 2y+N/2 ». U 2 +PV 



(3.37) 



Consequently (cf. Eq. p. 251) ). the sought radial momentum-space Coulomb Sturmian functions are 
found to be 



F nil (E,p) = 2"- l rli 



2 J y 7r(n r + 2Z + iV-2)! 



q {q 2 +p 2y+(N+l)/2 C n r W+p 2 ) 1 ' 

(it is seen that n r is simply the radial quantum number). They obey the weighted orthonormality 
relation 

/>oo 

JV— XfJi , „2\ 



dpp iV - V + q 2 )F nil (E,p)F n , l (E,p) = 5„ xK (3.39) 
and are standardized to be positive in the immediate vicinity of p = 0. 



3.5 Some properties of solutions to Eq. ( 13.51) 

Once the radial integral equation (|3.16|) is solved, we may return to the ./V-dimensional problem 
(J33J). 

It follows from Eq. (|3.36p that the degree of degeneracy of the eigenvalues in Eq. (|3.5p is higher 
than Sj N , as A„ r / depends on the quantum numbers n T and I through their sum. To distinguish 
between different eigenvalues, we introduce the principal quantum number 

n = n r + l + l, (3.40) 

in terms of which we have 

A„ r , = A n = . (3.41) 
Z <7b 

The degree of degeneracy of the particular eigenvalue A„ obviously is 

n-i 

D^ = J2^- 1} - (3-42) 

1=0 

In view of Eq. (|3.7p , this may be rewritten as 

1=0 v 7 z=o v 7 

The sums in Eq. f|3 .43[) are elementary; we have [37l Eq. (0.15.1)] 



fc=0 

and consequently 
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Comparison of Eqs. (|3.45[) and (|3.7[) shows that it holds that 

D^=dZ\. (3.46) 

(The relation in Eq. ()3.46j) is, of course, a symptom of the 0(N + 1) dynamical symmetry of the 
-/V-dimensional negative-energy Schrodinger-Coulomb problem.) 

Degenerate Sturmian eigenfunctions associated with the characteristic numbers (|3.41[) , labeled 
with the principal rather than the radial quantum number, are 

Z nlv (E,p) = F n ^ 1 j(E,p)Y l [ i N - 1 \n p ). (3.47) 

It is implied by Eqs. (|3.8[) and (I3.39[) that these functions obey the weighted orthonormality relation 

J d N p (p 2 + q 2 )Ki n {E,p)^n-v v -{E,p) = 6 nn ,6 w 6 vn >. (3.48) 

It remains to consider the completeness problem for the Sturmian set {E n ;^(iJ, p)}. It is known 

that the functions {(1 — f 2 )™/ 2- -^Cn" (£)}: with Rea > —4, form a complete set in the functional 
space L 2 ((— 1, 1)), the corresponding closure relation being 

2 2 - 1 [r(a)] 2 f> n r !(n r + a) {a) (a) 

ir ^ T(n r + 2a) K ^> l? ' [(l-e 2 )(l-e ,2 )]" /2 " 1/4 

(-1 < £, £' < 1, Rea > (3.49) 

with <J(£— £') denoting the Dirac delta distribution. We set in the above equation a = l + (N—l)/2, 
plug Eq. (|3.29p . combine the result with Eq. (|3.38j) . and then exploit the identity 

Hf-fl-S!^* ~P% (3-50) 

which follows from Eq. (|3.29l) and the well-known properties of the Dirac delta. This chain of 
operations transforms Eq. ()3.49j) into the closure relation for the radial Sturmians: 

V F n i(E,p)F n i(E,p / ) = 5<k P = ^ (3.51) 

Next, we have the closure relation 

d («-i> 

oo a i 

E E tf^'wC^K) = siN ' 1} ( n P - K)> ( 3 - 52 ) 

(=0 77=1 

which reflects the fact that the hyperspherical harmonics form a complete set in the functional 
space i 2 (S A ' _1 ). We insert Eqs. (|3.51[) and (|3.52[) into the right-hand side of the identity 

oo n— 1 oo oo 

EE E w^p)^(£,p') = E E ^Vp^r^K) E FnAE, P )F nxl (E,p') 

n=l 1=0 i7=l Z=0 17=1 n r =0 

(3.53) 

(3.54) 



and make use of the obvious fact that 



5 ] (P-P) = ( W /)(AT-l)/2 • 



This leads us to the symmetric closure relation 



-i<*r ij 



E E E w^p)^.*') = f" )( f v ^! if - ( 3 - 55 ) 

„=li=0 r,=l V(P" + 9 2 )(P 2 + 9 2 ) 

We have thus proved that the Sturmian functions {S„;^(S,p)} form a complete orthonormal 
set in the functional space L 2 2+q2 (M. N ); equivalently, the functions {\/p 2 + q 2 Yj n i v {E,p)} form a 
complete orthonormal set in the functional space L 2 (M. N ). 
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4 Determination of bound-state eigensolutions to the 
momentum-space Schrodinger— Coulomb wave equation 



d23]) 



We are now ready to determine bound-state eigensolutions to Eq. (|2.5p . following the procedure 
outlined at the beginning of Sec. 13.11 At first, we substitute A„, as given in Eq. (|3.41[) , for X a (E) 
in Eq. (|3. 21) . then use Eq. (13.41) . and solve the resulting equation for E. For each particular n, there 
is just one solution to that equation, which we shall denote as E n \ it is 



2 C 2 



E n 



Z 2 e 



2(n+^) (47re )a B 
Corresponding eigenfunctions, in accordance with Eqs. (|3.3I) and 

(n p ), 



(4.1) 



are 



r(N-l) 



with 

F n i{p) = AF n -i-ij(E n ,p). 

Equation (|2.5p shows that E~ Y is an eigenvalue of the Hermitian (in fact — real and symmetric) 
kernel 



(4.2) 
(4.3) 



2m^ (JV)(P ^ (4 7 re ll )2;T''^+ 1 ^/,. 



P 



,/|-(JV-l) 



From this and from Eqs. (|4.2I) and (|3.8p . one routinely deduces that the functions {$> n i v (p)} are 
orthogonal in the sense of 



/ d N p^ d Jp)$n'l'n'(p) ~ Snn'Sll'Srn'- 



(4.4) 



If, additionally, we require the functions under consideration are normalized to unity under the 
same scalar product that is used in the above equation, and if we choose the constant A in Eq. 
(14.31) to be real, we obtain the constraint 



(4.5) 



A 2 / d PP N ' 1 [F n ^ 1 . l (E n ,p)} 2 = 1. 
Jo 

To evaluate the integral in Eq. (|4.5|) , we use Eq. p. 381) . This yields 

N- 1 



d P p N ~ L [F n ^ hl {E n ,p)} 



■)2l+N-3 



r i 



(n-l-l)\(n+Z=Z) _ 2 
Tr(n + l + N-3)\ Qn 



x / i dMi-a z+A72 - i (i+6 l ) 



c 



(i+(iV-l)/2) 
n-l-l 



with 



q n = 



Zqt 



jV-3 
2 



(4.6) 



(4.7) 



It is evident that in the integrand on the right-hand side of Eq. (|4.6|) one may replace the factor 
1 + by 1 without altering the value of the integral. Applying then the formula (|3.33[) and 
imposing a further restriction on A to be positive, we obtain 



A = V2q n . 



(4.8) 



Hence, we conclude that the bound-state momentum-space eigenfunctions of the Schrodinger- 
Coulomb problem in K , N 2, associated with the energy levels (|4.ip and orthonormal in the 
sense of 

(4.9) 



/ d^P^nJ Jp)<$>n'l'T,'(P) = Snn'Sll>5 W , 
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may be chosen in the form 



2 / V 7r(n + ; + iV-3)! 

x n iV/2+l (4gnP)' r (i+(iV-l)/2) / gg-p 2 N y(^-l) r _ X / 4 1 

X9 « ((? 2 +p2) *+(iV + l)/2 »-*-l Ut+P 2 ; P 

The degeneracy of the energy level E n is, obviously, the same as that of the Sturmian eigenvalue 
A„, and is given by Eq. (|3.45|) . 

With Eqs. (|4.1[) and (|4.10l) being obtained, our task to determine bound-state eigensolutions 
to the Schrodinger-Coulomb problem in the TV-dimensional Euclidean momentum space is accom- 
plished. 
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Appendix: The Fourier transform of V(r) = Ar~ v ', < v < N, 
in R N 

In this Appendix, we shall follow closely Ref. [351 Sec. 5.9] (for an alternative approach, cf. Ref. 
Appendix 3]). We wish to evaluate the Fourier transform 



U(P) = * / d N re-^ h V(r) 



(A.l) 



where 

V(r) = — (0<v<N). (A.2) 



r 

To accomplish the task, at first we observe that 



/>oo />oc 



(A.3) 



hence, it follows that 

1 



r(|) 

If we plug Eq. (|A.4j) into Eq. (jA.lj) . we obtain 







U(n) = - f d N rc- [pr / h f°° d£ f/2-l e -«r 2 

(2^/^(1) A« y d44 



(27rft) w / 2 r (I) 

The integral over r is readily found to be 

N/2 



(A.5) 



/ d N rc-^/ fi e^ 2 = M e-P 2 / 4 ^, (A.6) 



so that Eq. (|A.5[) becomes 
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and further, after switching to the integration variable r\ = p 2 /4?i 2 £, 

(2h) N ^-» 1 f°° Hn JJV-,/)/2-l -, 

r(|) 



£f(p) = ^ ' / d^-^^e"". (A. 



n 



The remaining integral is elementary and we eventually arrive at 



In the case of the Coulomb potential 



Ze 2 

V ( r ) = -75-^. (A.10) 



(47re ) 



r 



from Eq. (|A~9|) we find that 



Ze 2 {2h) N ' 2 - 1 Y (^) 1 
(47re ) 



^(P) = -77TTT— ^ - 2 N-l - ( A ' n ) 
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